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approximation 



(N 
O 

o 



Thomas Bourdel * 
Laboratoire Charles Fabry, UMR 8501, Institut d'Optique, CNRS, 
Univ Paris Sud 11, Avenue Augustin Fresnel, 91127 PALAISEAU CEDEX 
(Dated: October 4, 2012) 

We study the superfluid transitions in bidimensional (2D) and tridimensional (3D) disordered 
and interacting Bose gases. We work in the hmit of long-range correlated disorder such that it can 
be treated in the local density approximation. We present the superfluid transition curves both 
in the disorder-temperature plane well as in the disorder-entropy plane in 2D and 3D Bose gases. 
Surprisingly, we find that a small amount of disorder is always favorable to the apparition of a 
superfluid. Our results offer a quantitative comparison with recent experiments in 2D disordered 
ultra-cold gases, for which no exact theory exists. 
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I. INTRODUCTION 

Disorder is an essential ingredient for the understand- 
ing of transport in condensed matter physics. At low 
temperature, it affects the conductivity of a metal and 
it even induces phase transitions to insulating states 
Superconductor-insulator transitions [ll are observed in 
thin metal films Q and, in high-Tc superconductors, 
where doping intrinsically introduces inhomogeneities in 
the CuO-planes Q. Understanding the complex inter- 
play between disorder and interactions in these systems 
remains a major challenge. Whereas the above men- 
tioned electronic systems are fermionic, superconductiv- 
ity originates from the bosonic nature of Cooper pairs. 
As long as disorder does not break the Cooper pairs, the 
problem is reduced to a study of dirty bosons . It has 
mainly been studied numerically in the framework of the 
disordered Bose- Hubbard model [H, 0| • This model was 
realized experimentally in ID and 3D ultra-cold gases in 
disordered lattices . The physics is comphcated by 

the presence of a Mott insulating transition, which occurs 
in the absence of disorder. 

In this paper, we focus on a different model with corre- 
lated disorder in a continuous Bose system. This model 
has been studied previously in 3D jn|. Our work is 
motivated by recent experiments studying the superfluid 
transition in disordered 2D ultra-cold atomic Bose gases 
[l^ [isj . A schematic phase diagram for such a system 
has been proposed and at low temperatures and low dis- 
order, a superfluid is expected [lj|. However, the exact 
shape of the superfluid region is not known. In this pa- 
per, we calculate the superfluid transition curve for a 
disorder that is correlated at long range such that we 
can make use of the local density approximation. In this 
limit, we find quantitative predictions in both 2D and 
3D. In 3D, we rely on a mean- field description of the 
Bose gas [la] whereas in 2D, we use the equation of state 
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from Monte-Carlo-simulations [15| . Interestingly, we find 
a very different quantitative behavior in 2D and 3D and 
that a small amount of disorder (either added reversibly 
or at constant temperature) is always favorable to the 
appearance of a superfluid. 

The paper is structured as follows. In the second sec- 
tion, we deflne our system parameters and the conditions 
of validity of the local density approximation. In the fol- 
lowing section, we calculate the critical disorder at zero 
temperature in homogeneous and harmonically trapped 
Bose gases in 3D and 2D. In the fourth section, the super- 
fluid transition curves are found as a function of disorder 
and temperature. In the fifth section, we extend our re- 
sults to get the curves as a function of entropy which is 
more convenient for the comparison to experiments which 
is presented in the sixth section. 



II. THE LOCAL DENSITY APPROXIMATION 
IN A DISORDERED POTENTIAL 



We consider Bose gases with repulsive contact interac- 
tion, i.e. with a positive scattering length a, such that 
there is no collapse of the wave function. The 2D gases 
are supposed to be harmonically confined in the verti- 
cal z direction with a frequency a;z/27r such that the 
atoms occupies in this direction a single quantum state, 
whose size is Iz = \J h/mujz, where h is the Planck con- 
stant and m the atomic mass. We also suppose that 
a <^ Iz, such that the interactions can be described 
by 3D scattering. The interaction parameter is then 
52D — {fi^/'m)g = {h^ /m)y/ST:a/lz 0|. In three dimen- 
sions, the interaction parameter is g^u = Ai:h^a/m. The 
healing length of a gas of density Tid, where d stands 
for the dimension, is ^ = y/mg^Y^n^. It reflects the 
size over which the density of a condensate (or quasi- 
condensate in 2D) comes back to its original value when 
a local perturbation is imposed. 

We consider a disorder created by a laser speckle, with 
a blue detuned frequency as compared to the atomic tran- 
sition frequency. In this case, the optical potential V{y) 
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is repulsive [l8|. We suppose that the speckle intensity 
has a gaussian correlation function whose 1 / radius is 
a. The mean disorder potential V is proportional to the 
laser power. Knowing a and V ^ all statistical properties 
of the disorder are then known JJi]. In particular, the 
probability density of having a disorder Y is an exponen- 
tial P{y) ~ exp(— y/y)/!^. The classical percolation 
threshold for such a disorder is proportional to the mean 
disorder V and we call a the proportionality coefficient. 
In 2D, a2D =0.52 and, in 3D, ago ~ 4 lO'^ [ll|. 

The much lower value of a in 3D as compared to 2D re- 
flects the fact that it is easier to find ways around the 
peaks of the disorder in 3D. 

The local density approximation (LDA) assumes that 
the gas can be considered locally as homogenous with a 
local chemical potential /^locai = A* ~ Kixt, where VJjxt is 
the external potential, i.e. the sum of the trapping and 
disorder potential. In typical experimental conditions in 
non-disordered harmonic traps, the LDA tends to give 
extremely accurate results. In order for this approxima- 
tion to be valid in a disorder potential, a must be much 
larger than all characteristic lengths of the gas such as 
the healing len gth ^ and th e thermal De Broglie wave- 
length AdB = \/ ^ttK^ / mk^T , where fee is the Boltzmann 
factor and T the temperature. 



III. SUPERFLUID TO INSULATOR QUANTUM 
PHASE TRANSITION 

In 2D and 3D, homogeneous interacting Bose gases at 
zero temperature are always superfluid independently of 
the density. Within the local density approximation, a 
superflow in the presence of disorder can then exist, if the 
region of non-zero density percolates through the sample. 
The quantum phase transition (T = 0) is thus reduced to 
the percolation of the density. The density of the gas 
can be found from the Thomas-Fermi approximation of 
the Gross-Pitaevskii equation [2^ which corresponds to 
the local density approximation for a (quasi-)condensate 
rid = Max(/iiocai, 0)/gdD- The gas will thus be superfluid 
when the region where /ziocai > percolates, i.e. when the 
global chemical potential ^ exceeds the classical percola- 
tion threshold of the disorder aV . For < /x < uzdV , 
local superfluid islands exist in the disorder minima but 
it does not lead to global superfluidity. 

In a disordered Bose gas at zero temperature in a box 
the number of atoms is : 



mean density 



N= / Max(^-l^(r),0)/.gdDdV 



(1) 



When averaging of the disorder, we can replace the posi- 
tion dependance of the disorder by an integral : 

{N)=[ [M(ixitx^V,0)/gdDd''rPiV)dV. (2) 
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. (3) 

which is valid in both 2D and 3D. The expansion for 
low value of ado is especially accurate in 3D. In 2D, we 
find the critical disorder strength V^i-it ~ 8.7g2D('^2D''^)j 
whereas in 3D, V^rit ~ 1-2 10''53d('^3d'^)- ^ much higher 
disorder is thus needed in order to reach the superfluid 
quantum phase transition in 3D. 

In the presence of an additional harmonic confinement, 
we suppose that the gas is superfluid as soon as the su- 
perfluid threshold is obtained in the center of the trap 
where the density is maximal. The calculation can then 
be done along the same lines with the trapping potential 
entering in the expression of the local chemical potential. 
The spacial integration leads to different results in 2D 
and 3D for the critical atom number as a function of the 
disorder : 

(N^D^) = vr^^(2 - 2cxp(-a2D) - 2a2D + ^L) 



(ATCondx 
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(4) 

(asD - w)^/^e-''dw 
(5) 

4V2 y5/2 2 7/2 



15 {hwfna^/^^^^LJJhJo 



Q3D^o 15 {hujf/'^a^Jmbj/h^ 



Here, lo is the geometrical mean of the two in-plane trap 
oscillation frequencies in 2D and of the three oscillation 
frequencies in 3D. 



IV. PHASE DIAGRAM OF 2D AND 3D BOSE 
GASES AS A FUNCTION OF TEMPERATURE 
AND DISORDER 

We now look for the temperature dependance of the su- 
perfluid transition in the presence of disorder. In the pre- 
vious section, we have found the critical disorder at zero 
temperature. The critical temperature with no disorder 
is known both in 2D [lB| and 3D [lB| ■ In the following, we 
study the shape of the superfluid transition curve which 
connects these two points in the temperature-disorder 
phase diagram. We separate the 3D case, in which a 
mean field approximation is possible, from the 2D case, 
in which one has to rely on classical field numerical cal- 
culations. 



The superfluid transition at finite temperature 
in a 3D gas 



The integration over the volume is then trivial. At the 
transition point, i.e. for = aV ^ we find the critical 



A 3D Bose gas at low temperatures is well described 
by two components, a condensate fraction and a thermal 
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fraction. For typical interaction parameters used in ex- 
periments, the interaction between the atoms in the ther- 
mal cloud can be neglected and the interaction between 
the thermal cloud and the condensate can be treated in a 
mean field approach. Similarly to the T = case, the dis- 
ordered 3D Bose gas is superfluid if the condensate den- 
sity percolates through the sample. The critical chemical 
potential is thus equal to a^uV. In order to know when 
the phase transition is reached, we have to calculate the 
relation between the atom number, the temperature and 
the chemical potential in the presence of disorder. 

The thermal cloud can be treated in the semi-classical 
approximation. The atom number for a Bose gas in a 
box is u&i 



th 



d^fcdV 



(27r)3 QP{h^ky2m+2g3onl'S^'^{r)+V{r)-f,) _ i 



(6) 

where /3 = l/fceT and where 2^30^3™'^ is the term due 
to the interaction with the condensate part in the mean 
field approximation. After averaging over the disorder, 
and integrating over r and fc, the mean thermal density 
is 



{nth)XlB = / e-^'PolyLog 



-, e 



where PolyLog [n,a] — ^ a'/r 



dv , (7) 
(8) 



1=0 



where the Polylog function is the result of the integration 
over k and where the absolute value comes from the mean 
field interaction with the condensate. 



The total mean density is then (ritot) 

^cond\ 



-,cond\ 



[nth), 

where {n'^^'^) comes from equation 3 and (nth) comes 
from equation 7. At the phase transition, the chemi- 
cal potential is /i = a^uV and the conservation of the 
density gives the relation between the disorder strength 
and the critical temperature. It can be written in a di- 
mensionless form using the quantities V' = V / (ntot)53D, 
r = {{ntot)XlBr'^', and a' = {{ntot)a'y/^ which 
quantifies the interaction strength. The equation relating 
V' and T' at the phase transition is then 



1 = V'{{~1 + exp(-a3D) + aso)) 



+T' 



3/2 



e '"PolyLog 



-2\a3u~v\a'V' /T' 



(9) 

dv . (10) 



It can be solved numerically. 

The boundary between the normal and superfluid 
phase for a disordered Bose gas in a box is plotted in fig. 
[TJ For low disorder, we find that the critical temperature 
increases with disorder. This is due to the compression 
of the gas by the disorder, which leads to an increase 
of the effective phase space density. When the disorder 
approaches its critical value at T = corresponding to 
V' « 1.210^, the increase of the percolation threshold 
starts to play a role and the critical temperature quickly 
goes to zero. 




0.01 



T' 



FIG. 1: Superfluid phase transition of a disordered 3D Bose 
gas in a box. We have fixed {ntot)a'^ ~ 10~®. At low temper- 
ature and low disorder the gas is superfluid. 



A similar calculation, albeit slightly more technical, 
can be done in the presence of an harmonic confine- 
ment. In this case, natural non-dimensional param- 
eters are T' = kBT/huN^/^, V = V/fujN^/^, and 
the strength of the interaction is characterized by a' = 
N^^^ay^muj/h.The equation of the superfluid boundary 
is then 

l = r{a,^~vf'^e-^dv 
15a' Jo 



/ ( y^-e-''Erfi[V^^) 
Jo 



[3 V 
xPolyLog -,exp(--|a3D - wl — ) 

where Erfi is the imaginary error function. 



dv 



(11) 




FIG. 2: Superfluid phase transition of a disordered 3D Bo se 
gas in an harmonic trap. We have fixed a = N^^^a^moj/h = 
0.023. At low temperature and low disorder the gas is super- 
fluid. 

The boundary between the normal and superfluid 
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phase for a disordered Bose gas in an harmonic trap is 
plotted in fig. [5J The interpretation is the same as be- 
fore except that the disorder induced enhancement for 
the critical temperature is not as strong because heating 
the gas results in a further spreading in the trap and in 
a reduced density. 



B. The superfluid transition at finite temperature 
in a 2D gas 

The physics of 2D gases is more complex than the 
one of 3D gases because the effects of interactions are 
much stronger and cannot be captured in a mean field 
approximation. The superfluid phase transition at flnite 
temperature is of Berezinskii-Kosterlitz-Thouless type 
[ISHjISEl and is intrinsically a beyond mean- field effect. 
As a consequence it is not possible to separate the gas into 
two components. The phase space density of the homo- 
geneous 2D Bose gas D = 7120 A^g has to be a function of 
fj,/T and g because the interaction parameter g is dimen- 
sionless. This equation of state is not analytically known. 
It has been computed in Monte-Carlo simulations close 
to the phase transition . Further from the transition, 
simple analytic forms are available. At low phase-space 
density the mean-field approximation works, and at high 
phase space density the Thomas-Fermi limit is rapidly 
approached. Numerically, it is possible to interpolate be- 
tween the different regimes with a good accuracy. 

From Monte-Carlo studies, we know that a finite su- 
perfluid fraction is expected at a critical value of the ratio 
Mcrit/^B?^ ~ qLog(13.2/ff)/7r which depends on the inter- 
action parameter g |26j . In the local density approxima- 
tion, we expect the gas to be superfluid if the regions of 
flnite superfluid fraction percolate, i.e. if the difference 
of chemical potential /i — /icrit is larger than aV. To flnd 
the phase transition for a Bose gas in a box, the equation 
to solve is thus : 



"2d)A^b= P{V)D{(3{fi,,it+a2BV~V),g)dV . (12) 



Using the non dimensional parameters V' — 
V / g2T){n2T)) , T' = l/(n2D>AdB, and g, the phase dia- 
gram takes the form shown in fig. [31 At ^ = 0, the phase 
transition takes place at (n2D)A^B « 8 ~ 1/0.125 as pre- 
dicted in the absence of disorder for g — 0.96 [l^. A 
T = 0, we recover the moderate critical value expected 
from the treatment using the Thomas-fermi approxima- 
tion. In between, the shape of the phase diagram is differ- 
ent from the 3D case. We interpret this as a consequence 
of the much higher percolation threshold in 2D. In the 
phase diagram in 2D, there is no region where the perco- 
lation threshold is irrelevant in contrast to the 3D case. 
Nevertheless at small disorder, the effect of compression 
is still dominant and the critical temperature increases 
with disorder albeit by a small amount. 

Again, a similar calculation can be made in an har- 
monic trap. In this case, we chose V = V/hujN^/^ and 
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FIG. 3: Superfluid phase transition of a disordered 2D Bose 
gas in a box. g — 0.096. At low temperature and low disorder 
the gas is superfluid. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
T' 

FIG. 4: Superfluid phase transition of a disordered 2D Bose 
gas in a harmonic trap, g = 0.096. At low temperature and 
low disorder the gas is superfluid. 



T' = ksT/hajN^/'^ in addition to g as non-dimensional 
variables. The phase diagram is then shown in fig. 21 The 
initial increase of the critical temperature with disorder 
is still present although it is a relatively weak effect. 



V. CRITICAL ENTROPY AS A FUNCTION OF 
DISORDER 

In all the previous cases, the addition of a small 
amount of disorder is favorable to the apparition of a 
superfiuid fraction, in the sense that the critical temper- 
ature increases. This is due to the fact that the disor- 
der compresses the gas in a smaller region and therefore 
leads locally to an increase of the phase space density 
(at constant temperature). In a recent experiment with 
ultra-cold atoms |12| , the disorder is added reversibly af- 
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ter the preparation of the gas. In this case, the entropy 
is conserved (as weh as the number of particles) but the 
temperature increases due to the effective compression of 
the gas. In the following, we derive the phase diagram as 
a function of disorder and entropy per particle, which is 
more directly related to the experiment. 

In the local density approximation, the calculation of 
the critical entropy can be done as soon as the entropy 
of a homogeneous gas is known. In the 3D be- 
fore the semi-classical approximation can be used and the 
entropy per unit of volume is [l^ 



^ = ^BAd^(-PolyLog 



:i pWocal/fcBT 



/^local 



PolyLog 



— e' 
2' 



,Mlocal/fcBT 



(13) 



<S'2D = fcB(2Prcd/-D — / k^T) , whcrc the reduced pres- 
sure Prcd = J-!^^^ Dd{fi/kBT) can be calculated from 
the equation of state. 





S[kB] 

FIG. 5: Superfluid phase transition of a disordered 3D Bose 
gas in a box as a function of entropy per particle and disorder. 
(fitot)a^ = 10~^. At low entropy and low disorder the gas is 
superfluid. 

This expression can be used to find the phase diagram 
as a function of the entropy and the disorder. It is plot- 
ted in fig. [5] Without disorder, the critical entropy per 
particle is the one that is found at the Bose-Einstein con- 
densation threshold in a box. When increasing the dis- 
order, the critical entropy first increases to a value close 
to 2.1. This is a consequence of the change of the shape 
of the density of state at low energy [27] from propor- 
tional to \/E without disorder to proportional to E"^^^ 
with disorder. Only at a disorder of the order of the crit- 
ical disorder at T = 0, the critical entropy decreases. In 
a 3D experiment, although it may appear counterintu- 
itive, it is thus possible to induce superfluidity in slowly 
adding some disorder. A similar statement holds in the 
3D harmonically trapped Bose gas, although we do not 
show the plot in this paper. 

In the 2D case, the calculation of the entropy is pos- 
sible as soon as the equation of state is known be- 
cause of the scale invariance in 2D interacting Bose 
gases [2^. More precisely, the entropy per particle is 
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FIG. 6: Superfluid phase transition of a disordered 2D Bose 
gas in a harmonic trap as a function of entropy per particle 
and disorder, g = 0.096. At low entropy and low disorder the 
gas is superfluid. 

In fig. ini the phase diagram of a harmonically trapped 
2D gas is plotted as a function of entropy per particle 
and disorder. For small disorder, the critical entropy 
still increases with disorder and it is thus favorable for 
the appearance of a superfluid fraction. However, this 
effect due to the change in the density of state is rapidly 
hindered by the effect of the increase of the percolation 
threshold in contrast to the 3D case. 



VI. COMPARISON TO EXPERIMENTS 

Experimentally, the appearance of coherence have been 
studied in ultra-cold 2D rubidium disordered gases p^ . 
In this experiment, we expect noticeable deviations from 
the LDA because cr ^ w AdB- Nevertheless, since there 
is no quantitative theory in the regime of microscopic 
disorder explored in the experiment, it is interesting to 
compare with the LDA approximation results. The two 
disorder strengths used in the experiment correspond to 
V « 0.3 and V « 0.7. The value V = 0.7 is close to 
the critical disorder and the transition is strongly shifted. 
This is in agreement with the observation of no sudden 
appearance of phase coherence in the range of entropy 
explored in the experiment. The value V = 0.3 corre- 
sponds to a lowering of the critical entropy per particle 
by 0.3 /cb- This value is somewhat larger than the shift 
^ 0.2 /cb found in the experiment for the sudden appear- 
ance of coherence, a phenomenon linked to the superfluid 
phase transition. Note that the absolute values of the en- 
tropy per particles differ between the above theory and 
the experiment because of a small fraction of atoms that 
populates the excited states of the vertical strongly con- 
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fining harmonic oscillator. We have checked that adding 
the vertical levels does not significantly affect the value 
of the entropy shift in the LDA approximation. 

Another quantity that we can calculate in the local 
density approximation is the change of temperature due 
to the addition of the disorder. For the parameters of 
the experiment, the temperature changes are found to be 
5nK and 14 nK. These values clearly exceed the exper- 
imentally observed heating. The fact that the LDA ap- 
proximation overestimates the effects of the disorder can 
be interpreted as a consequence of an effective smoothing 
of the disorder for a ~ ^ ~ XdB and thus as a proof of 
beyond LDA physics. 

VII. CONCLUSION 

In this paper, we have analyzed 2D and 3D Bose gases 
in the presence of disorder, which is treated in the local 
density approximation. The superfluid part of the phase 
diagram is found as a function of disorder and temper- 
ature and as a function of disorder and entropy. Our 



results show the competition between two effects when 
adding disorder. One is the appearance of a percolation 
threshold, which tends to hinder superfluidity. The sec- 
ond is the disorder induced compression of the gas and 
the change of the density of state. On the contrary, it is 
favorable for the appearance of superfluidity. Since the 
percolation thresholds greatly differ in 2D and 3D, we 
find very different quantitative behaviors of the super- 
fluid transition curves. Although it might be counterin- 
tuitive, the addition of a low amount of disorder always 
lowers the critical temperature as well as the critical en- 
tropy. Finally, this works permit quantitative compar- 
isons with experimental findings in the presence of disor- 
der. 
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